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Atmospheric escape from HD 189733b observed in Hi Lyman-a: 
detailed analysis of HST/STIS September 2011 observations 
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ABSTRACT 

Observations of transits of the hot giant exoplanet HD 189733b in the unresolved Hi Lyman-Q line show signs of hydrogen 
escaping the upper atmosphere of the planet. New resolved Lyman-a observations obtained with the STIS spectrograph onboard 
the Hubble Space Telescope in April 2010 and September 2011 confirmed that the planet is evaporating, and furthermore 
discovered significant temporal variations in the physical conditions of its evaporating atmosphere. Here we present a detailed 
analysis of the September 2011 observations of HD 189733b, when an atmospheric signature was detected. We present specific 
methods to find and characterize this absorption signature of escaping hydrogen in the Lyman-a line, and to calculate its 
false-positive probability, found to be 3.6%. Taking advantage of the spectral resolution and high sensitivity of the STIS 
spectrograph, we also present new results on temporal and spectro-temporal variability of this absorption feature. We also 
report the observation of HD 189733b in other lines (Sim at 1206.5 A, Nv at 1240 A). Variations in these lines could be 
explained either by early occultation by a bow-shock rich in highly ionized species, or by stellar variations. 
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1. Introduction and observations 

1.1. Evaporation of hot Jupiters 

Although th e phenomenon of evapora tion had been an- 
ticipa ted by Burrows fc Lunin^ ( 1995 ) and Guillot et al 
(jl996[) . the existence of a large number of hot Jupiters 
surviving atmospheric escape suggested that evapora- 
tion processes should be modest. In this frame it 
came as a surprise to dis cover that the exoplane t 



HD 209458b was losing gas (jVidal-Madiar et all 120031 ). 



Transit observations of the Lyman-a line showed ex- 
cess absorption due to an extended cloud of neutral hy- 
drogen Hi and constrained t he (H i) escape rat e with 



a lower limit of 10 



10 



(Vidal-Madi ar et al.l 12003 : 



Vidal-Madiar fc Lecavelier des Etangsij2004f). T his discov 



ery has been challenged by iBen-Jaffell ([20071), but the 
apparent discrepancy has been resolved and the result 
obtained by Ben-Jaffel on thi s first data set strength- 
ens the evaporation scenario (jVidal-Madiar et al.l 12004 ; 
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Vidal-Madiar et al.l 120081) . It was also confirmed by two 



subsequent observations at low spectral resolution: a 5% 
absorption of the whole Lyman-a line was indeed mea- 
sured using unresolved emission line fiux during plane- 
tary transits observed with the Spac e Telescope Imaging 
Spect rograph (STIS) instrument ( Vidal-Madiar et al.l 



20041 ) and the Advanced Camera for Surveys (ACS) in- 

20081) onboard the Hubble 



strument ( Ehrenreich et al 



Space Telescope (HST). An i ndependent analysis of the 
low-resolution data set used bv lVidal-Madiar et al.l (|2004[ ) 
has confirmed that the transit depth in Lyman-a is sig- 
nificantly greate r than the transit depth due to th e plane- 
tary disk alone ( Ben-Jaffel &: Sona Hosseini 2010l ). These 
three independent observations show a significant amount 
of gas at velocities exceeding the planet escape velocity, 
leading to the conclusion that HD 209458b is evaporat- 
ing. This conclusion is strengthened by the observations 



of ab sorption in the O i and C ii lines (jVidal-Madiar et al 



20041 ) , which led to the identification of a blow-off escape 



mechanism. This is also confirmed by recent observations 



2 



V. Bourrier et al.: Atmospheric escape observed in Hi Lyman-a 



of absorption in the C ii and Sim lines using th e HST 



Cosmic Origins Spectrograph ( Linskv et al. 201Cll ). Note 



that additional analyses of these observations are being 
made (Ballcster & Ben-Jaffel, in prep.). A similar con- 
clusion has been reached in the case of WASP- 12b with 
observati ons of Mg il and nea r UV broad-band transit ab- 



sorption (iFossati et al. 

Based on these observational constraints, several 
theoretical models have been developed to explain and 
characterize the evaporation pro c esses (iLammer et al. 



Baraffe et al.l 
Jaritz et al.l 



2004 . 
20051: 



2005 



2003t iLecavelier des Etangs et all |2004 l2007l. l2008b 



2007: Holmstrom et al. 



Tian et al.l 
20081: 



et al. 

Murrav-Clav et al.l 12003: lAdamsl l201lt IGuoI I2OIII) . In 1-2. HD 189733b 



m)A _ Ye\\4 2004, 2006' 



sit observations of HD 189733b made in April 2010 and 
September 2011 with the UV channel of the repaired 
HST/STIS (using the sam e methodology as employed by 
Vidal-Madiar et all [2OO3I for HD 209458b), and showed 
that the escape from this planet is subject to significant 
temporal variations. In this paper we present a detailed 
analysis of the STIS transit observations of HD 189733b 
made in September 2011. We also show a new timc- 
resolved spectral signature of the absorption observed at 
this epoch by using time-tagged photon count data. We 
also present the transit observed in other lines, such as 
Sim at 1206.5 A and Nv at 1239 A and 1243 A. 



the models dev e loped by lHolmstr5m et al.l (120081 ) 



20051: 

Stone fc Progal 



Garcia Munoz 
2009t 



Ekenback et al 



In 

and 
pro- 



( 2010f ). the observed hydrogen is 
duced through charge exchange between the stellar wind 
and the planetary escaping exosphere. Even in these last 
models, neutral hydrogen escaping the planet is required 
for charge exchange with hot protons of the stellar wind. 
Therefore, all models led to the conclusion that most of 
the EUV and X-ray energy input by the host star is used 
by the at mosphere to escape the planetary gravitational 
potential ( Owen fc Jacksonll2012; for detailed estimates o f 



the EUV and X-ray input, see ISanz-Forcada et aLll2011l) . 
With this idea in mind, ILecavelier Des Etangsl (|2007l) 
developed an energy diagram in which the potential 
energy of the planet is plotted versus the stellar energy 
received by its upper atmosphere. For HD 209458b, 
the escape rates derived from this diagram are not 
high enough for the evaporating exoplanet to lose a 
signif icant amount of its mass (jEhrenreich fc Desert 
l201l[ ). However, the nature of planets with smaller 
masses and closer orbits could be significantly altered 
by evaporation, and these planets would end as plane- 

The 
and 



2004) 



20091) 



tary remnants (jLecavelier des Etangs et al . 
recently disc overed CoRoT - 7b ([Leeer et al 
Kepler- 10b ( Batalha et al.l I2OIII) could be examples 
of a new category of planetary remnants, which were 
proposed for classification as " chthonian planets" (see 
Lecavelier des Etangs et al. I l2004l) . 



For a few years, the only observed evaporating planet 
was HD 209458b, and many questions remained to be an- 
swered. What is the evaporation state of other hot Jupiters 
and very hot Jupiters? Is evaporation a common process 
amongst hot Jupiters? How does the planetary system and 
stellar characteristics influence the escape rate? Some light 
has been shed o n these questions th anks to the discovery 
of HD 189733b (|Bouchv et al.ll2005l ). a planet transiting 
a bright and nearby K star (V=7.7). Among the stars 
harboring transiting planets, HD 189733b belongs indeed 
to the second brightest at Lyman-a. Despite the failure 
of the HST/STIS in strument one year before the dis- 



(l20ln^ 



covery of HD 189733b. lLecavelier Des Etangs et al. 
detected atmospheric escape through Lyman-a transit 
observations with the HST / ACS instrument . Moreover , 
Lecavelier des Etangs et all ( 20121) compared two tran- 



Thc very hot Jupiter HD 189733b is located 19.3 parsecs 
away from Earth, with a semi-major axis of 0.03 AU and 
an orbital period of 2.2 days. The planetary transits and 
eclipses ca n be used to probe its atmosp here and environ - 



ment {e.g., Charbonneau et al. 2008, : Desert et al. 2009t ). 



HD 189733b orbits a bright main-sequence star, with a 
magnitude V=7.7. This K2V star emits one of the high- 
est Lyman-a flux ever measured for tra nsits observations 
(with the noticeable exception of 55 Cnc: lEhrenreich et al 



20121 ). which makes HD 189733b a particularly good can- 
didate to study atmospheric evaporation. 

Thanks to the proximity of HD 189733, and the 
large surface covered by the planet against its rela- 
tively small star, extensive observational studies have 
been made of this hot Jupiter. The planet has a 
mass Mp=1.13 Jupiter masses (Mjup) a nd a radius 



i?P= 1 .16 Jupiter radii ( R.^uu) in the visible (|Bakos et al 



20061 : 

suits m a 



Winn et al.l l2007l ). The planetary disk transit re- 



2.4% occultatio r i depth from vis ible to near 
infrared (|Desert et al.ll2009t ISing et al.llioilh . The short 
period of the planet (2.21858 days) has bee n mea- 



sured precisely (iHebrard fc Lecavelier Des Etangsl 12006 



Knutson et al.l l2009f ). Spectropolarimetry has measured 
the strength and topolog y of the stellar magnetic field 



which reaches up to 40 G (jMoutou et al ] |2007t iFares et al, 
2OIOI) . Sodium has been dete cted in the pla.net a tmo- 



sphere by both ground-b ased (iRedfield et al.l l2008l ) and 



space-base d observatio ns ( Huitson et al.ll2012[ ). Using the 
HST /ACS. IPont et al.l (|2008l ) detected atmospheric haze, 
whic h is interpreted as Mie scattering by small parti- 
cles (ILecavelier Des Etangs et al.ll2008al ) . High signal-to- 



noise HST/NICMOS/ST IS observations dsing et al.l l200d 
201 ll ) and HST/WFC3 (jCibson et aDl2012l ) haveshown 
that the near-IR spectrum below 2 / xm is due to haze 



scatterin g and/or wat e r abs orption (jSwain et al.l 12008 



see also iGibson et al.l 120111) . and Rayleigh scatterin; 



migh t even extend to longer wavelengths (|Pont et al.l 
20121). The dete ction of an H2O s ignature using tran- 



sit photometry ( Tinetti et al.l 20(37) h a s been subjected 
to controversv (lEhrenreich et al.l 120071: lAgol et al-lboiOl: 
Desert et al.ll2009il201ll) ." It has been tentatively proposed 
that CO molecules can explain the excess absorption seen 
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at 4.5 ^tm ( Charbonneau et al.l[2008 : Desert et al. 2009t 
Knutson et al.1 12012) . Applying best-estimate approaches 
to dayside infr ared emission spectr a of HD 189733b, 
Lee et all (|2012[ ) and lLine et~all (|2012l) both reported the 
detection of H2O and CO2 in the atmosphere (see also 
Swain et al. 20091). Using Spitzer spectroscopy of plane- 
tary eclipses. lGrillmair et al.l ( 2008 ) also found evidence of 



H2O absorption signatures and possibly weath er-like vari- 
ations in the atmospheric conditions. Note that lAgol et al.l 
( 2010( ) found an upper limit of 2.7% on the variability 
of the dayside planet flux, which rules out the most ex- 
treme weather fluctuations on HD 189733b. An alternative 
source for the observed variability may be the change of 
location of a daysid e hot spot detected at an offset froin 
the substellar point jKnutson et al.ll2007ilAgol et al.ll2010l 
Maieau et al.l 20121 ). A sensitive search with GMRT has 
provided very low upper limits on the meter-wavelength 
radio emission from t he planet, indicating a weak plan- 



etary magnetic field ( Lecavelier Des Etangs et al. 20091 



20111). 



i 



1.3. Previous Lyman-a observations of HD 189733b 



2. Observations 

The observations of HD 189733b in the Lyman-a line per- 
formed in 2007 - 2008 showed significant absorption in 
the exosphere of the planet, while in April 2010 no atmo- 
spheric absorption signature was detected. To address the 
question of these variations, new transit observations of 
the planet were made in September 2011 with the G140M 
grating of the HST/STIS. The log of these time-tagged 
observations is given in Table |T] The analysis described 
in this paper was made with ID spectra extracted us- 
ing CALSTIS (version 3.32) data pipeline. These spectra 
cover the far-UV wavelengths from 1195 to 1248 A with a 
spectral resolution of about 20kms~^ at 1215.67 A. The 
planet was observed during four consecutive orbits of the 
HST around the Earth: two orbits before the transit, one 
during the transit, and one after the transit. The first ex- 
posure lasts 1800 s and the three others about 2100 s. Data 
acquisitions are interrupted by Earth occultation lasting 
about 3500 s. Contrary to ACS observations, for which the 
co-addition of three independent transits was necessary 
to detect the transit signature in the Lyman-a line, the 
highest sensitivity of the STIS data allows an independent 
study of each transit observation. 



In 2007-2008, because of the failure of the Space 
Telescope Imaging Spectrograph (HST/STIS) in 2004, 
HD 189733 was observed with the Solar Blind Camera of 
the (HST/ ACS). Transit observations showed a diminu- 
tion of 5.05ib0.75% of the entir e Hi Lyman-a curve 
( Lecavelier Des Etangs et al. 2010t ). This was more than 
the 2.4% planetary disk occultation depth and was inter- 
preted as the result of the evaporation of the planet's up- 
per atmosphere. However, the limited spectral resolution 
of the ACS data in the far-ultraviolet wavelengths of the 
Lyman-a line (at 1215.67 A) implies that no conclusion 
could be made on the radial velocities of the escaping gas, 
and this detection called for more observations. After the 
refurbishment of HST in May 2009, HD 189733b was thus 
observed during transit in April 2010 with the G140M 
grating of the repaired HST/STIS. Disappointingly, no 
significant abso rption was detected in the Lyman -a line 
(2.9±1.4%, see iLecavelier des Etangs et al.l |2012|) . other 
than the o ccultation de p th by the planetary disk alone. 
Note that IJensen et al.l (|2012[ ) reported a significant ab- 
sorption feature at Ha in the transmission spectra of 
HD 189733b. This transit-dependent absorption was de- 
tected within a narrow band at the line center (-24.1 
to 26.6 km s~^). Because of the limited resolution of the 
ACS spectra, Jensen et al. could not directly compare 
the velocity ranges between their H-a absorption detec- 
tion and the H i Lyman-a absorpti on feature reported by 
ILecavelier Des Etangs et"aL ( 2010t ). Although the stellar 
Lyman-a emission line is strongly absorbed by the inter- 
stellar hydrogen over their observed velocity range, our re- 
cent and highly resolved STIS observations could enable 
a more direct comparison between the two atmospheric 
signatures. 



3. Data analysis 

3.1. Resulting spectra 

We identified in the STIS spectra the stellar emission 
lines of Sim (1206.5 A), Ov (1218.3 A), the Nv doublet 
(1242.8 A and 1238.8 A), and the bright Hi Lyman-a line 
(1215.67 A) (Fig. [T|). To calculate the radial velocities rel- 
ative to the star, we need to estimate the radial velocity 
of the star relative to the STIS wavelength calibration. 
Using the stellar emission lines we found that the star has 
a red-shifted velocity of ■~5kms~^ (heliocentric). 

In order to detect possible transit signatures of 
HD 189733b and its atmosphere, we calculated for each 
stellar line the time evolution of the total emission flux, 
or of the flux within a given wavelength range. We also 
compared the April 2010 and September 2011 observa- 
tions of the stellar lines other than the Lyman-a line 
(for comparison of this line betw een 2010 and 2011, see 
Lecavelier des Etangs et al.ll2012[ ). A signature of the at- 
mosphere is considered to be detected when an excess ab- 
sorption is seen in addition to the absorption by the plan- 
etary disk itself during the planet transit. 

3.2. A/v doublet (at 1238.8 A) and Ov line 

For the April 2010 and September 2011 observations, we 
compared the total flux in the stellar Sim emission line 
and the Nv doublet before, during, and after the transit. 
Results are summarized in Table We do not detect any 
variation showing a transit-like signal in the 2010 or in the 
2011 observations for the Nv line at 1238.8 A. All varia- 
tions, in particular differences in flux between 2010 and 
2011, are most likely attributable to stellar variations in a 
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Data set 


Date 


Observation 


Time from center of transit 






Start 


End 


Start 


End 


HST orbit #1 


2011-09-07 


19:49:21 


20:20:55 


-03:31:02 


-02:59:28 


HST orbit #2 


2011-09-07 


21:18:37 


21:56:10 


-02:01:46 


-01:24:13 


HST orbit #3 


2011-09-07 


22:54:26 


23:31:59 


-00:25:57 


00:11:35 


HST orbit #4 


2011-09-08 


00:30:16 


01:07:49 


01:09:52 


01:47:25 



Table 1. Log of the September 2011 observations. Time is given in UT. 
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Fig. 1. Plot of the wavelength-cahbrated STIS spectrum 
of HD 189733b with the Lyman-a, Si ill, O V, and N v stel- 
lar emission lines. 

line arising from an active region of the stellar atmosphere. 
The analysis of the O V line, fainter and located within the 
Lyman-a line red wing, gave no significant results. 

3.3. Sim line and A/v doublet (at 1242.8 A) 

Although there is no significant variation in the first epoch 
observations for the Sim line and the Nv doublet, in the 
second epoch we see significant time variations of the flux 
measured in the Sim line, and in the Nv line at 1242.8 A 
(see Table [5]). In the 2011 observations, these stellar lines 
are observed to be fainter before the transit than during 
and after the transit (Fig [2] and Fig [3]). There is a 3.5% 
probability to find the high variations observed between 
the fluxes before, during, and after the transit for the Nv 
line at 1242.8 A (see the values in Table [5]). This prob- 
ability is 0.6% for the Sim line. Variations in the Sim 
line brightness have alread y been observed in the plane- 
tary system of HD 209458b (jLinskv et al.ll2010h . However, 
here we see a higher flux during the transit, which can- 
not be interpreted by a classical absorption in the plan- 
etary exosphere. Our observation could be interpreted 
by an early ingress caused by the formation of a bow- 
shock surro unding the planet magr ietosphere. Th i s was 



p t „ 

modeled bv lLlama etldl (l201lh and lVidotto et all (|2010l . 
2011allbh . and explains for example the WASP- 12b obser- 



vations of an early ingress in the Mg ii line (jFossati et al 



2010l Has well et al. 2012[) . Alternatively, intrinsic stellar 



line variations could be the explanation of our observa- 
tions. This raises the question of possible stellar line vari- 
ations triggered by the planet, as already suspected in the 
case of the C ii flare detected in ACS observat ions of the 
same planet ( Lecavelier Pes Etangs et al.ll2010[) . The Sim 
line and Nv doublet thus appear as possible candidates 
to diagnose the star- planet interaction (SPI), and addi- 
tional observations are needed to address this possibility. 
Hereafter we consider only the Lyman-a observations. 




1206.0 



1206.5 
Wavelength (A) 



1207.0 1206.0 



1206.5 
Wavelength (A) 



1207.0 



Fig. 2. Plot of the flux in the Sim line for the 2010 
(left panel) and 2011 (right panel) observations. The black 
dashed line shows the spectrum before the planet transit; 
the green solid line with error bars at the 1 a level shows 
the flux during the transit; the red dotted line shows the 
flux after the transit. 



3.4. Lyman-a line 

3.4.1. The stellar Lyman-a emission line profile 

The Lyman-a line is the brightest stellar line in the STIS 
spectra from 1195 to 1248 A. The single stellar Lyman-a 
line is massively absorbed by the interstellar atomic hy- 
drogen Hi in a narrow band at the line center (1215. 6A), 
and to a lesser extent by the interstellar deuterium D I 
around 1215.3 A. As a result, the spectrally resolved line, 
at the resolution of the G140M spectrograph, is composed 
of two peaks separated by this deep absorption. In order 
to estimate the Lyman-a line flux and profile, we fitted 
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Sim 


inc 




N V doublet 






Flux within 1206 - 1207 A 


Flux within 1238.4 - 1239.1 A 


Flux within 1242.5 - 1243.1 A 




(10"^*ergs 


-icm-2) 


(10^^^ erg 


s ^ cm ^) 


(10-^5 


erg s ^ cm ^ ) 




2010 


2011 


2010 


2011 


2010 


2011 


Before transit 


1.05±0.04 


0.87±0.04 


3.59±0.19 


2.77±0.16 


1.74±0.13 


0.97±0.10 


During transit 


0.94±0.05 


0.98±0.05 


3.15±0.23 


2.72±0.22 


1.50±0.16 


1.40±0.15 


After transit 


1.07±0.06 


1.08±0.06 


3.07±0.23 


2.77±0.22 


1.75±0.17 


1.19±0.14 




3.82 


10.1 


3.86 


0.04 


1.65 


6.70 



Table 2. Total flux in the Sim emission line (at 1206.5 A) and the N v emission line (at 1242.8 A and 1238.8 A) before, 
during, and after the transit. The values are the standard deviations associated to each stellar line in either 2010 
or 2011. 




1242.5 1243.0 
Wavelength (A) 



1242.5 1243.0 
Wavelength (A) 



Fig. 3. Plot of the flux in the N v hue between 1242.5 and 
1243.1 A, for the 2010 (left panel) and 2011 (right panel) 
observations. The black dashed line shows the spectrum 
before the planet transit; the green solid line with error 
bars at the 1 <j level shows the flux during the transit; the 
red dotted line shows the flux after the transit. 



the line profile as observed i n Sep temb er 2011, following 
the m ethod of Wood et al. ( 2005 ) and lEhrenreich et al. 
(|201l[) . The best result was obtained by modeling the 



stellar emission line using two Voigt profiles with the 
same width and damping constant. The free parameters 
of the fit are the total stellar fiux, the Voigt parameters 
(width, damping constant, and the wavelength difference 
between the two profiles centers), and the ISM parameters 
(Hi column density, ISM and star radial velocities, and 
t urbu lence velocity). We used a D/H ratio of 1 .5x10"^ 
(e.g.. iHebrard fc MoodbOOl iLinskv et al.ll2006l) . Finally 
we added a single free parameter to fit the strength of the 
wing of the line spread function (LSFV which is added to 
the G140M published tabulated LSItl- The wing of this 
LSF is taken as a Gaussian with a fixed width measured 
in the pubhshed G140M LSF to be (JwingLSF — 4.6 pixels. 
The comparison between the best model reconstructed 
Lyman-a stellar profile and the observed spectrum yields 
a of 36.1 for 49 data points in the wavelength range 
1214.25-1215. 5A and 1215.8-1217.1 A (see Fig. S]). The 



Model 




k 


DOF 


BIG 


AIG 


Reference model 


36.1 


9 


40 


71.1 


54.1 


TiSM 7^ 


35.8 


10 


39 


74.7 


55.8 


Double Gaussian LSF 


34.4 


12 


37 


81.1 


58.4 


2 different Voigt profiles 


35.5 


11 


38 


78.3 


57.5 


1- Voigt profile 


41.5 


8 


41 


72.6 


57.5 


2-Gaussian profile 


166.6 


8 


41 


197.7 


182.6 



Table 3. x^j BIG and AIG for various models with k 
degrees of freedom used to fit the Lyman-a line profile. 



resulting profile shows a double peak emission, which is 
usually seen in the Lyman-a line of cool stars, and can be 
inter preted as self-abso rption in the chromosphere of the 
star (jWood et al1l2005[) . 

In the search for the best model to fit the Lyman-a 
profile, we used the Bayesian Information Griterion (BIG) 
and the Akaike In formation Griterion (AIG) (see, e.g., 
de Wit et al. l l2012l) . These criteria prevent us from over- 



fitting and are based on the likelihood function given by 
the x^ s-nd on ^ penalty term rel ated to the number of 



^ Ihttp: //www. stsci . edu/hst/stis/perf ormance/spectral. 



parameters in the fitting model (jGrossfield et al.l 12012 : 
Gowan et al ]l2012h . The lowest BIG is obtained for the ref- 
erence model described above (see parameters in Table [3]). 
When we allowed for the temperature of the ISM to be 
a free parameter (allowing different turbulent widths for 
the ISM hydrogen and deuterium absorption lines), we 
obtained a better x^ but the larger BIG and AIG show 
that this decrease cannot be interpreted in terms of the 
information content of the data. We also obtained a bet- 
ter fit to the data by using a two-Gaussian LSF with three 
additional free parameters, but again the larger BIG and 
AIG shows that the increase of freedom is responsible for 
the apparent improvement in the fit. We also concluded 
that the better fit to the profile obtained by using two dif- 
ferent Voigt profiles (with different widths and damping 
constants) is not significant. We also tested if more simple 
emission line profiles with lower degrees of freedom could 
be used (one single- Voigt profile or two-Gaussian profile) 
and concluded that the reference model with two similar 
ri^aigtt profiles provides the best fit to the data. 
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Fig. 4. Plot of the theoretical profile of HD 189733 
Lyman-a line. The black thin line shows the theoretical 
intrinsic stellar emission line profile as seen by hydrogen 
atoms escaping the planetary atmosphere. The black thick 
line shows the resulting profile after absorption by the in- 
terstellar hydrogen (f2f5.6A) and deuterium (1215. SA). 
The hue profile convolved with the HST G140M instru- 
mental line spread function (red line) is compared to the 
observations (blue histogram), yielding a good fit with a 
of 36.1 for 40 degrees of freedom. 



3.4.2. Geo-coronal emission 

In the raw data, the stellar emission line is superimposed 
with the geo-coronal air glow emission from the up per at- 
mosphere of the Earth (jvidal-Madiar et al.ll2003l) . This 
emission line can be well estimated and removed from the 
final spectrum using the calstis data pipeline (version 2.32 
of November 5, 2010). Independent re- analysis of raw data 
using the s ame method o logy as in IVidal-Madiar et al. 
(|2003l) and lOesert et aP t004) confirmed that the air- 
glow emission can be efficiently subtracted and included 
in the final error budget. Moreover, because we used a 
narrow spectrograph slit of 0.1", the airglow contamina- 
tion is limited to the central part of the Lyman-a line and 
does not contaminate the part of the spectrum where at- 
mospheric signatures are detected, in the wings of the line 
(Sect. E]). The airglow emission may be affected by daily 
and seasonal variations, and the data of September 2011 
present a noticeably low airglow emission level (Fig. [S]). 



3.4.3. Telescope "breathing" 

Using time-tagged acquisitions of the STIS spectra, we 
checked variations at short time scal es within a given H ST 
orbit in a similar manner as done bv lBen-Jaffef(l2007l) . In 
the 2011 observations, we see a periodic variation of the 
measured Lyman-a flux in phase with the HST position 
on its orbit, with the same amplitude of about 10% at all 
wavelengths and no spectral signature. These variations 
are not seen at a significant level in the 2010 observations. 
They arc interpreted by the STScI STIS experts as caused 
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Fig. 5. Plot of HD 189733b Lyman-a line profile (black 
line) superimposed with the geo-coronal emission (red and 
blue lines). The red line corresponds to the 2010 observa- 
tions, the blue line to the 2011 observations. All lines are 
averaged over the four HST orbits. The geo-coronal emis- 
sion is noticeably low during the 2011 observations. 



by the telescope "breathing", which is known to affect 
STIS observations in the UV and in the visible because 
of variations of the telescope throughput. The breathing 
is caused by changes in temperature that the HST ex- 
periences during each orbit. This causes small motions 
of the secondary mirror in the STIS bench , which alters 



the a lignment of the target in the aperture (iKimble et al 
19981 ). It may also change the focus position and the size 



of the point spread function, affecting the spectrograph 
throughput through the narrow 0.1" slit. Although the 
breathing depends on a large number of parameters and 
cannot be predicted, it can be modeled f or a past ob- 
servation. Using the STScI Focus model (jCox fc Niemi 



20111 ). we thus assessed the repeatability of the breath- 
ing effect on the focus position during the 2011 obser- 
vations. While the time windows of these observations 
correspond to significant variations of the focus position, 
the 2010 observations coincide with times when the fo- 
cus position is subjected to only low variations (Fig. [5]). 
The shape and amplitude of these systematic variations 
can change between visits of the same target, as was 
found in STIS opti c al transit observations of HP 209458 
( Brown et al. 2001 : Charbonneau et al] 20021: Sing et al 



2008 ) and of HD 189733 (|Sing et al.ll2011l: iHuitson et al 



20121) ■ and in STIS FU V transit observations of 55 Cnc 
([Ehrenreich et al~ 2012 ). Extensive experience with opti- 
cal STIS data over the last decade furthermore indicates 
that the orbit-to-orbit variations within a single visit are 
both stable and highly repeatable, exce pt for the first or- 
bit, which displays a different trend (see Brown et al.ll2001 



and ICharbonneau et al.l 120021 ). Because of its repeatabil- 
ity, this effect is easy to correct for. 

For each HST orbit, we divided the time-tagged expo- 
su res into six 379 s sube xposures. Following the method 
of ISing et al.l (|201l[ ) and I Brown et al.l (|200l[ ). we phase- 
folded these exposures over the HST-orbital period of 
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Fig. 6. Longitudinal motion of the secondary mirror for 
the 2010 observations (top panel) and the 2011 observa- 
tions (bottom panel) . Time is given relative to the center 
of the planetary transit. Hatched zones show the time win- 
dows of our observations of HD 189733b. 



96 minutes. Each exposure was then spectrally summed 
over the entire Lyman-a line, and normalized by the mean 
value of all the exposures of its related orbit. We excluded 
the first orbit. Wc also excluded observations obtained 
during the planetary transit. Wc found that the breath- 
ing effect was best fitted with a se cond-order polynomial 
(Fig. [TJ see lEhrenreich et al. 2012 for a similar correction 
of this effect). 




0.1 0.2 
HST orbital phase 

Fig. 7. Lyman-a fluxes for 379 s exposures, phase-folded 
on the HST orbital period. The empty circles correspond 
to the first orbit. The green circles correspond to the 
measurements obtained during the planetary transit. The 
black line is the second-order polynomial fit to the other 
orbits measurements (filled black dots). 

Fig.[7]shows that the first orbit seems to be lower than 
the expected unabsorbed Lyman-a stellar light curve. 
This could mean the absorption signatures described in 
Sect, m and obtained with the two orbits before the tran- 
sit as reference, may be underestimated. We correct the 
data for the breathing effect by dividing the spectra of 
each 379 s subexposure by the value of the polynomial 



fit at the same HST orbital phase. This correction does 
not apply to the first orbit because of its different trend. 
The spectra are reconstructed for each orbit by coadding 
their respective exposures (Fig. [S]) . Note that we obtain 
the same estimates for absorptions and ratios measured 
using the full HST o rbit exposures, with or without th e 
correction, as done bv lLecavelier des Etangs et al.l (j2012l ). 
However, the correction must be taken into account for 
measurements on partial orbit exposures obtained using 
time-tagged data. 



1.15 
1.10 
I 1-05 

T3 

0) 

:^ 1.00 
CO 

E 

I 0.95 
0.90 
0.85 



-4-3-2-10 1 2 3 

Time (hours) 

Fig. 8. Raw fluxes not corrected for the breathing effect 
(empty squares), as a function of time in the September 
2011 observations. Time is given relative to the center of 
the planetary transit. Circles correspond to the corrected 
flux. Each point stands for a 379 s exposure spectrum, 
summed over the whole Lyman-a line, and normalized 
by the flux over the entire corresponding orbit. Vertical 
dashed lines show the beginning and end of ingress and 
egress of the transit. 



4. Atmospheric hydrogen 

4.1. Full Lyman-a flux 

In a flrst step, to search for the transit signature of the 
planetary atmosphere in the Lyman-a line, we compared 
the total line flux measured during the transit with the 
flux measured before the transit. The transit depth by the 
planetary disk alone as seen from the visible t o the near- 



infrared amounts to 2.4% ( Pont et al.l l2007t ISing et al 



20091 : iDesert et al.ll2009t ). The observations of HD 189733b 
made in April 2010 showed no significant excess in the 
transit absorption depth, with a depth for the flux within 
the whole Lyman-a li ne of 2.9zbl.4%, including the plan - 
etary disk occultation (jLecavelier des Etangs et al.ll2012l ). 
In contrast, a transit absorption of 5.0±1.3% is found 
in the total flux of the Lyman-a line for the 2011 ob- 
servations. We fitted the fiux measured during the tran- 
sit with a classical planetary occultation model from 
Mandel k Agoll (|2002[) . The best fit is consistent with the 
presence of an extended absorbing atmosphere surround- 
ing the planet (see Fig. [51). Substracting the planetary disk 
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absorption, the signature measured during the 2011 tran- 
sit corresponds to an excess absorption, due to atmosperic 
hydrogen only, of 2.3±1.4%. This is consistent with the 
excess absorption of 2.71±0.75% due to hydrogen only, 
obtained with the HST/ACS obs ervations in 2007-2008 
( Lecavelier Pes Etangs et al. 2010h . 




Fig. 9. Plot of the total flux within the whole Lyman-a 
line, as a function of time (blue square symbols). Each 
HST orbit is divided into two 1127 s independent subex- 
posures, corrected for the breathing effect (except for the 
first orbit, see Section [3.4.3p . Time is given relative to the 
center of the planetary transit. Vertical dashed lines show 
the beginning and end of ingress and egress of the tran- 
sit. Horizontal error bars centered on the symbols show 
the duration of the exposures. A light curve fitted with 
a classical planetary occultation model is displayed as a 
solid black line, with a transit depth in excess compared to 
the transit depth observed at optical wavelengths (dashed 
black line). The bottom panel shows the corresponding 
residuals (empty symbols). 



4.2. Spectrally resolved absorption features 

In contrast to the ACS spectra, the spectrally resolved 
STIS spectra allow us to look for absorption signatures 
within specific wavelength intervals (Fig. [TO)) . Two clear 
absorption regions can be seen within the Lyman-a line in 
the 2011 transit observation. The most significant absorp- 
tion signature is visible in the blue part of the line from 
-230 kms"^ to -140kms~^. Another feature, although less 
significant, is visible at the peak of the red wing from 
BOkms""'^ to 110kms~^ (Fig.[TT|). There seems to be some 
absorption within the same velocity region in the blue 
wing for the post-transit observation, albeit with a lower 
depth. The red region displays no absorption feature in 
the post-transit observation. In the following sections we 
describe the methods used to quantitatively characterize 
these two signatures, in particular their precise velocity 
range and absorption depths. 
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Fig. 10. Lyman-a line profile. Four spectra, correspond- 
ing to the four HST orbits of the September 2011 observa- 
tions, are displayed as a function of radial velocity relative 
to the star. The black and blue lines show the fiuxes be- 
fore the planet transit; the green line with error bars at 
the 1 (T level shows the flux during the transit; the red line 
shows the flux after the transit. 
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Fig. 11. Zoom on the blue wing (left panel) and red wing 
(right panel) of the Lyman-a line profile displayed in 
Fig. [To] The black and blue dashed lines show the fluxes 
before the planet transit; the green solid line with error 
bars at the 1 a level shows the flux during the transit; the 
red dotted line shows the flux after the transit. 



4.3. Quantitative estimates 

To quantitatively characterize the absorption features 
within the Lyman-a line, we studied the ratio of the flux 
during transit to the flux before transit. The transit sig- 
nature is the combination of the absorption by the at- 
mospheric neutral hydrogen and the 2.4% occultation of 
the stellar disk by the planetary disk at all wavelengths. 
To characterize the absorption due to the hydrogen only, 
we decreased the flux before the transit by a factor cor- 
responding to the planet absorption. Then, we searched 
for the strongest signature at every possible wavelength 
range, excluding intervals narrower than 0.1 A (about two 
pixels in the STIS spectra). We looked for the most sig- 
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nificant absorption, characterized by the highest value of 
the ratio of the relative absorption to its noise (Fig. [T2|) . 

In the 2011 transit observations, the most significant 
signature is found within the range -230 to -140 km s^^, 
yielding an Hi absorption depth of 12.3±3.7% (3.3 cr de- 
tection). This signature, due to the atmospheric hydrogen 
only, is the most significant over the entire Lyman-a pro- 
file. The measured velocity of the gas is higher than the 
escape velocity of 60kms^^ from the planet (the effec- 
tive escape velocity is even lower than GOkms"-'^ at higher 
altitude where we observe the absorption). The total ab- 
sorption depth of 14.4±3.6%, which includes the planetary 
disk occultation, also corresponds to high altitudes: it is 
equivalent to a disk with a radius of 2.8 Jupiter radii. 
Therefore, the observed neutral hydrogen must be at high 
altitudes, escaping the planet gravity. Note that if we ex- 
clude the first HST orbit in the calculation of the fiux be- 
fore the transit, the absorption depth is the same within 
la (18.3±5.6%), but the detected range is somewhat nar- 
rower, from -230 to -180kms~^. 
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Fig. 12. Plot of the fiux (upper panel), absorption (bot- 
tom panel), and ratio of the absorption to its noise (middle 
panel) in the blue part of the Lyman-a line, as a function 
of velocity. The upper panel displays the fiux before the 
transit multiplied by the occultation factor of the plane- 
tary disk (black dashed line), to be compared to the flux 
during the transit (blue solid line). In the lower and mid- 
dle panels, total cumulative values are calculated within 
domains of increasing range. The red dotted line is for ra- 
tios calculated within domains between -230kms~^ and 
increasing velocities. The green dash-dotted line is for ra- 
tios between -140kms~^ and decreasing velocities. 



The same method was applied to the red wing of the 
Lyman-a line, with a search window in the range 40 to 



200 km s^^. We flnd the most significant transit absorp- 
tion signature in the range 60 to 110 km s"^ (see Fig. [13]). 
Substracting the occultation by the planetary disk, this 
signature yields a 5.5±2.7% absorption depth {2.0 a de- 
tection) . This detection is the same if we exclude the first 
orbit, whether the spectra are corrected for the breathing 
effect or not. Cumulative transit depths due to the plan- 
etary disk and the atmospheric hydrogen can be found in 
Table H 




100 150 
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Fig. 13. Same as in Fig. [T2]in the red part of the Lyman- 
a line. The red dotted line is for ratios calculated over 
domains between 60kms~^ and increasing velocities. The 
green dash-dotted line is for ratios between 110 km s~-^ and 
decreasing velocities. 



We also divided each exposure of a given HST orbit, 
corrected for the breathing effect, into two independent 
spectra with 1127 s exposure time. We measured the ab- 
sorption over the range -230 to -140kms~^ for the April 
2010 and September 2011 observations. A light curve fit- 
ted to the data using a model of a single planetary disk 
occultation unambiguously shows that, in contrast to the 
2010 measurements, the 2011 flux measurements require 
the presence of an extended hydrogen atmosphere. Each 
of the two independent transit measurements of 2011 arc 
consistent with a 14.4±3.6% absorption, interpreted as at- 
mospheric absorption in addition to the planetary disk oc- 
cultation fFig.rn]). Note that the absorbing atmosphere is 
likely not spherical, which explains why the post-transit 
data of the 2011 observations cannot be fitted with the 
simple model used here (see Section |476)) . A full descrip- 
tion of the geometry of the hydrogen cloud is beyond the 
scope of the present paper, however. 
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Transit depths 


HST/ACS (2007 - 2008) 
not spectrally resolved 


HST/STIS (2010) 
spectrally resolved 


HST/STIS (2011) 
spectrally resolved 


Planetary disk occultation 


2.4% 


2.4% 


2.4% 


Whole Lyman-a line 


5.1±0.8% 


2.9±1.4% 


5.0±1.3% 


-230 to -140 km s"^ 




0.5±3.8% 


14.4±3.6% 


60 to llOkms"^ 




5.8±2.6% 


7.7±2.7% 



Table 4. Transit depths in a given wavelength range (including planetary disk occultation and atmospheric absorp- 
tion), for the three observations of HD 189733b. 
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Fig. 14. Plot of the flux between -230 and -140kms~^ in the blue wing of the Lyman-a line as a function of time relative 
to the center of the planetary transit (filled symbols). The red triangular symbols are for the 2010 observations, while 
the blue square symbols correspond to the observations of 2011. Each HST orbit is divided into two 1127 s independent 
subexposures, corrected for the breathing effect (except for the first orbit, see Section 13.4. 3p . Vertical dashed lines 
show the beginning and end of ingress and egress of the transit. Horizontal error bars centered on the symbols show 
the duration of the exposures. The light curve of the planet transit at optical wavelengths is displayed as a dashed 
black line. A light curve fitted with a classical planetary occultation model to the flux during the transit is displayed 
as a solid black line, and shows that the excess absorption feature detected in the 2011 data is not seen in the 2010 
data. The bottom panel shows the corresponding residuals (empty symbols). 



4.4. False-positive probabilities 

To assess the significance of the absorption features found 
in the Lyman-a light curve, we made the hypothesis that 
the transit spectrum shows no absorption and calculated 
with a bootstrap-like method the false-positive probability 
to have signatures as significant as the one we detected, 
but caused by noise only. 

From a reference-averaged spectrum, we generated two 
random spectra simulating spectra before and during the 
transit free of atmospheric absorption. The two spectra 
were obtained by adding random Gaussian noise to the 
reference spectrum, with an amplitude corresponding to 
the error estimated for the spectra before and during the 
transit. We generated these two spectra for 1,000,000 runs, 
and in each run used the absorption detection method de- 
tailed in Section H31 to find the most significant absorption 
signature. Due to a better consideration of the searched 
velocity intervals, we found l ower estimates of the false- 



positiv e probability than in iLecavelier des Etangs et al 



( 2012h . Indeed, taking into account that the range of 
the searched-for absorption feature cannot be narrower 
than the spectral resolution of the spectrograph, we ex- 



cluded false-positive signatures found on wavelength inter- 
vals covering 1 or 2 pixels, which would be narrower than 
the STIS resolution. The false-positive probability to find 
an excess absorption depth at more than 3.3 a in the range 
-350 to -50kms~^ is only 3.6%. It is unlikely that this par- 
ticular signature comes from statistical noise in the data. 
The false-positive probability to find an excess absorption 
depth at more than 2.0 cr in the range 40 to 200kms~^ 
is 24.6%. With this higher probability, the 5.5±2.7% ab- 
sorption depth detected in the red wing is possibly due 
to noise in the data. Note, however, that a similar signa- 
ture was also observed in the case of HD 209458b, with 
an absorption depth of 5.2zfcl .0% in the red win g be- 



tween 1215.89 and 1216 .43 A (jVidal-Madiar et al. 
Vidal-Madiar et al.l[2008l) 



2003 



4.5. Spectral and temporal variability 

4.5.1. Absorption depth as a function of radial velocity 

Here we address the possibility to resolve the absorp- 
tion feature in the blue range of the Lyman-a line profile 
between -230 and -140 kms"^. Although the absorption 
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depth seems to be higher for higher negative velocities, 
a Hnear fit to the spectral absorption within this domain 
shows that the slope is not significantly different from 0, 
with a value of (-1.5±1.3) x lO'^ ( km/s)-i (FiglHl). We 
conclude that we do not see significant structures in the 
profile of the H i absorption. 
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Fig. 15. Plot of the spectral absorption in the blue wing 
of the Lyman-a line in the 2011 observations (blue his- 
togram). The black line is the best linear fit to the ab- 
sorption and is not significantly different from a constant 
value. 



4.5.2. Variations of absorption with time during the 
transit 

The spectral resolution and high sensitivity of the STIS 
data allow the search for possible differences in the velocity 
structure of the absorption as a function of time during the 
transit. Indeed, the escaping hydrogen ma y form a curved 
cometary tail trailing behin d the planet (jSchneiter et al. 
2007tlEhrenreich et al.l2008[) . Since the tail crosses in front 
of the star, we can use time-tagged observations to probe 
different areas of the atmospheric cloud by measuring the 
absorption profile at different phases of the planet transit. 
We divided each exposure of a given HST orbit, corrected 
for the breathing effect, into three spectra with 751 s ex- 
posure time. Wc measured the absorption over the range 
-230 to -140 kms^^. The absorption over low negative ve- 
locities (-175 to -140kms~^) seems to decrease with time 
during the transit, while we see an increase of absorption 
for higher negative velocities (-230 to -175kms~^) (see 
Fig.lHl). 

For the first of the three exposures (shortly after the 
end of the ingress) the difference between absorption at 
low velocities (21±8.1%) and high velocities (9.6±8.9%) 
is not significant. However, if we coadd the last two expo- 
sures to obtain the absorption at the center of the transit, 
we find that the absorption at low velocities is lower than 
for the first exposure, with a value of 6.6±6.6%, while it 
has increased at high velocities, with a value of 24±5.9% 
(Fig. [T7|) . This is a 2tT difference between the absorption 
depth at high and low velocities. Numerical simulations 
of the dynamics of the atmospheric cloud crossing in front 
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Fig. 16. Plot of the total absorption over two complemen- 
tary velocity intervals, as a function of time during the 
transit, for three spectra with 751 s exposure time. Green 
disks stand for -230 to -175 km s^^ and black circles for 
-175 to -140 kms^^. Horizontal bars indicate the duration 
of each exposure. 
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Fig. 17. Plot of the total absorption over two complemen- 
tary velocity intervals, as a function of time during the 
transit, for two spectra with respectively 751 s and 1502 s 
exposure times. Green disks stand for -230 to -175kms~^ 
and black circles for -175 to -140kms~^. Horizontal bars 
indicate the duration of each exposure. 



of the star may help explain these variations (Bourrier et 
al. in prep). More observations of HD 189733b transit will 
also be necessary to confirm possible variations of absorp- 
tion with velocity over time. 

4.6. Post-transit observation 

Since the extended planetary atmosphere may form a 
curved cometary tail trailing behind the planet, the transit 
of the escaped hydrogen can last longer than the transit of 
the planet itself. Observations obtained after the transit 
may thus be affected by the signature of escaping atmo- 
spheric gas in the line of sight toward the star. This is 
the reason why, in the present analysis, we compared the 
Lyman-a line profile observed during the transit to profiles 
observed before the transit only, excluding post-transit ob- 
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servations. The best velocity interval -230 to -140kiiis~^, 
found for the in-transit observations, still yields an absorp- 
tion depth of 6.83±3.82% in the post-transit light curve 
(Fig nil). 

Following the method detailed in Sect. 14. 3[ we looked 
for the velocity interval providing the most significant ab- 
sorption in the post-transit spectrum. Wc found the most 
significant absorption within the blue wing of the Lyman- 
a line in the range -230 to -180kms~^, with a signature 
of 9.8±5.4%. The higher velocities, with regard to the ve- 
locities of the transit signature detected in the transit 
spectrum, might be interpreted as absorption by atoms 
accelerated to their maximum velocity at the end of the 
comctary tail. 
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Fig. 18. Plot of the flux measured within the range -230 
to -140 km s^^, for three 751 s subexposures per orbit dur- 
ing the transit and after the transit, as a function of time 
(blue square symbols) . Each flux, corrected for the breath- 
ing effect, is normalized by the flux of the corresponding 
exposure of the second orbit (i.e., the orbit just before 
the transit). Vertical dashed lines show the beginning and 
end of ingress and egress of the transit. Horizontal error 
bars centered on the symbols show the duration of the ex- 
posures in each HST orbit. The light curve of the planet 
transit at optical wavelengths is displayed as a solid black 
line. 



5. Conclusions 

Wc observed the hot Jupiter HD 189733b with the 
HST/STIS instrument in the far-UV wavelengths from 
1195 to 1248 A in April 2010 and September 2011. While 
no significant variation was detected in the Nv line at 
1238.8 A for both observations, the second epoch observa- 
tion shows absorption in the Sim line and the Nv dou- 
blet at 1242.8 A, albeit before the transit. This cannot be 
interpreted by a classical absorption in the planetary ex- 
osphere, and may be due either to an early occulation by 
a bow-shock or to intrinsic stellar line variations. 

We also presented a detailed analysis of the Lyman-a 
line of neutral hydrogen Hi observed in September 2011. 
A strong atmospheric absorption feature of 12.3±3.7% (in 



addition to the 2.4% absorption by the planetary disk) 
was found between -230 to -140kms~^. The false-positive 
probability to find such a signature in the blue wing of 
the line was revisite d to be only 3.6%. Since H D 189733 is 
an active star (e.g., ISanz-Forcada et al.ll201ll) . one might 
wonder if stellar activity might have been the source of this 
absorption signature. However, this feature is observed 
only within a very specific wavelength range, during the 
planet transit. The velocities of the absorbing hydrogen 
atoms, greater than the highest estimate of the escape ve- 
locity from the planet (60kms~^) and explained by the 
interaction between stellar wind and hydrogen escaping 
the planet, unambiguously show that neutral hydrogen 
was escaping from HD 189733b in 2011. 

HD 189733b has been observed in the Lyman-a line 
during five transits at various epochs. The three tran- 
sit spectra observed with the HST /ACS ins trument in 
2007/2008 (|Lecavelier Des Etangs et all |2010[) displayed 
time variations in the transit depth at ~2.5 a. The absorp- 
tion depth measured in the HST/STIS spectra observed 
in September 2011 is consistent with the absorption depth 
in one of the three ACS spectra, and confirms the hydro- 
gen escape from the planetary atmosphere. On the other 
hand, no absorption was found in the HST/STIS spectra 
observed in April 2010, which is consistent with the last of 
the three ACS spectra. Roughly half of HD 189733b tran- 
sit observations thus show detectable escape processes in 
neutral hydrogen. 

Here we also report possible differences in the velocity 
structure of the absorption as a function of time during 
the transit. The Lyman-a line also shows some absorption 
in the range -220 to -180 km s^^ in the post-transit spec- 
trum. This is consistent with the idea of a cometary tail 
of accelerated hydrogen atoms trailing behind the planet. 
Upcoming observations of HD 189733b in the Lyman-a 
line, along with a model of the hydrogen escape and its ac- 
celeration mechanisms, will help understand the temporal 
and spectral variability of the absorption in the exosphere 
of the planet. 
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